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We present results on the sum of the masses of light and strange quark using improved staggered quarks. Our 
calculation uses 2+1 flavours of dynamical quarks. The effects of the dynamical quarks are clearly visible. 



1. INTRODUCTION 

1.1. Motivation 

Simulations of hardronic matter with realistic 
sea quark effects form a long standing problem 
in lattice gauge theory. To achieve this goal on 
the upcoming generation of computers, improved 
staggered quarks are most promising. Compared 
to fermion formulations based on the Wilson ac- 
tion, staggered quarks have fewer variables and 
are faster to simulate for small quark masses. 

Naive staggered quarks suffer from large flavour 
changing strong interactions, which can be sup- 
pressed by the use of fat links ffl . This action can 
be further improved to ©(agO , a**) 0, the most 
precise to date, see 1^ for recent developments. 
These actions have proven to be successful to re- 
duce the splittings in the pion spectrum [Q and 
the large renormalisations |^,^, as observed for 
naive staggered quarks. Due to the advantages of 
improved staggered quarks it has been possible 
to generate data sets with 2-1-1 flavours of dy- 
namical sea quarks with a low TO^/mp-ratio 
An alternative proposal for the improvement of 
staggered quarks is the hypercubic action 

Staggered quarks have a residual "chiral" sym- 
metry, which protects against additive mass 
renormalisations. The point of vanishing quark 
mass is exactly known beforehand. This removes 
the dependency on chiral extrapolations in the 
determination of the strange quark mass. 



1.2. Simulation data 

For our calculation, we use the simulation re- 
sults of the MILC collaboration j?). They use 
a fermionic action corrected to 0(01^0^ ^a^) 
which is referred to as ASQTAD. There are re- 
sults for 0, 2 and 2-1-1 flavours of dynamical 
fermions. Their lightest quark mass gives a low 
mT^/nip « 0.37. The simulation parameters of the 
conflgurations have been tuned to keep the lattice 
spacing fixed in units of the gluonic observable ri . 

2. LATTICE SPACING 

2.1. Upsilon spectrum 

T-spectroscopy on the configurations of is 
discussed in ||l^ . Eor the data set with n/ = 2-1-1 
and the smallest m^/mp ratio one gets a^^ ~ 
1.6 GeV. In the quenched approximation the in- 
verse lattice spacing is substantially larger. 



2.2. Kaon spectrum 

In real world, the difference my — "^ps be- 
tween the square of the mass of the vector and the 
pseudo-scalar meson is almost independent from 
the masses of the valence quarks. In principle this 
should yield a good method for determining the 
lattice spacing. Vector mesons are unstable how- 
ever. Their decays will eventually be seen in dy- 
namical lattice results, but currently it is unclear 
what effect this has on the measured masses. 
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than the real world.. For data sets with dy- 
namical fermions on the other hand the split- 
tings rrip — mj, m|-. — m|- and m| — m^^ dif- 
fer by only 2%. On the Uf = 2 + 1 sets one 
gets — 1.4 GeV from mj^, — m\. For the 
quenched data this procedure gives an inverse lat- 
tice spacing which is about 100 MeV smaller. 

2.3. Gluonic scale 

The gluonic observable ri is determined in a 
similar way to rg, but gives the point where 
r^-F'(r) = 1. Since the configurations have been 
generated for fixed ri, it is natural to use ri as 
an alternative to set the scale. MILC estimated 
ri = 0.35 fm, using rg = 0.5 fm. Neither of 
these quantities is easily related to a physical ob- 
servable, however, and so there must be consider- 
able uncertainty in their physical value. Taking 
ri = 0.35 fm gives a^^ — 1.5 GeV for all configu- 
rations. 
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Figure 1. Interpolation of the pseudoscalar mass 
to the Kaon 



3. QUARK MASS 

3.1. Determination of the bare quark mass 

We determine the sum of the bare masses of 
the light and strange quark by interpolation of 
the pseudo-scalar meson mass to the physical K 
with the ansatz 



(amps)^ ~ a + b arrii^y + c(ar7iav)^ 



(1) 



With TOav := ^("^1 + 'n^2) we denote the average 
of the valence masses. This is shown in figure ^ 
for the most chiral dynamical data set. The data 
points represent the tt, K and rjs- The difference 
between a quadratic and linear interpolation is 
just visible in the upper right corner of the figure. 
At the physical K this difference is completely 
irrelevant. 

At this point the advantage of not having ad- 
ditive mass renormalisations for staggered quarks 
pays back. No simulation is needed to determine 
the point of vanishing quark mass. 

3.2. Renormalisation and matching 

The renormalisation and matching to the MS- 
scheme is done in a three step procedure. We use 
the quark self-energy to convert from the lattice 
mass mo to the pole mass |^ . In a second step we 
convert to the MS-mass at the scale of the pole 



mass and finally run to the scale /i = 2 GeV. At 
1-loop level we get 
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(2) 



("2") 

Here b := drriQ /mo — 2/tt log(amo). This is the 
1-loop contribution to the mass shift between the 
bare lattice mass and the pole mass with the loga- 
rithmic divergence taken out. For small values of 
amo this becomes a constant . For n/ = 2 + 1 
we obtain 



0.12 < Z^^^ < -0.04, 



(3) 



depending on the quantity used to fix the lattice 
spacing. It is interesting to note, for nai've stag- 
gered quarks, that is without fat links to suppress 
the flavour changes, Zq « 3. 

For as we use the average of ap for the scales 
aq = 1 and 2, determined on these configs [ill] . 
The difference between the two is included into 
the error. Preliminary calculations of the modi- 
fied BLM scale aq* |12] give values lying between 

(2) 

1 and 2. Because of the small Zq , variations in 
q have only very small effects. 
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Figure 2. Sum of light and strange quark mass 
in the MS-scheme at a scale fi — 2 GeV vs the 
r7i7r/mp-ratio of the lightest sea quarks. 



Table 1 

Preliminary results on m\ 



■ m, in MS-scheme. 



a{K*'^ ~ K^) a(T(lP) - T(15)) 




2 

2+1 



94(9) MeV 
83(11) MeV 
78(14) MeV 



71(7) MeV 
71(10) MeV 
71(13) MeV 



ical light quarks and one flavour of dynamical 
strange. This has been possible due to the use 
of improved staggered quarks. Due to this im- 
provement, the renormalisations are small and 
well controlled. The dynamical quarks clearly re- 
duce the dependency on the quantity used to de- 
termine the lattice spacing. Compared to other 
lattice calculations and sum rule calculations, see 
Q resp. Q for review, our results favour a low 
value of the strange quark mass. 

Acknowledgement 

This work is supported by the EU IHP pro- 
gramme, NSF, NSERC and PPARC. 



Our final result for the sum of the masses of the 
light and the strange quark in the MS-scheme is 
presented in figure ^. Different shapes of the sym- 
bols refer to the numbers of sea quarks. The dif- 
ferent shadings of the symbols indicate the deter- 
mination of the lattice spacing. The dashed line 
corresponds to the real world m7r/mp-ratio. The 
error bars encompass statistical uncertainties, the 
uncertainties of the determination of a and the 
uncertainty of the scale setting of as- One needs 
to add an uncertainty for the 2-loop effects. At 
the present time we estimate this as 0{a1) « 20% 
for Uf = 2 + 1 and smaller otherwise. 

In the quenched case, the results show a wide 
spread, depending on how a is determined. This 
spread is clearly reduced for Uf — 2 and even 
smaller for Uf = 2 + 1. There might even be a 
trend for a decrease with decreasing m-^/mp. In 
table we give preliminary results from the most 
chiral point of the corresponding data set. 

4. SUMMARY 

We presented the first ever calculation of the 
light quark masses using two flavours of dynam- 
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